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Abstract 

The Skyrme model can be generalised to a situation where static fields are maps from one 
Riemannian manifold to another. Here we study a Skyrme model where physical space is two- 
dimensional euclidean space and the target space is the two-sphere with its standard metric. The 
model has topological soliton solutions which are exponentially localised. We describe a superpo- 
sition procedure for solitons in our model and derive an expression for the interaction potential 
of two solitons which only involves the solitons' asymptotic fields. If the solitons have topological 
degree 1 or 2 there are simple formulae for their interaction potentials which we use to prove the 
existence of solitons of higher degree. We explicitly compute the fields and energy distributions for 
solitons of degrees between one and six and discuss their geometrical shapes and binding energies. 



1 A Skyrme Model in Two Dimensions 



The Skyrme model is a non-linear theory for SU{2) valued fields in 3 (spatial) dimensions 
which has soliton solutions. Each soliton has an associated integer topological charge or 
degree which Skyrme identified with the baryon number A soliton with topological 
charge one is called a Skyrmion; suitably quantised it is a model for a physical nucleon. 
Solitons of higher topological charge, called multisolitons, are classical models for higher 
nuclei. 
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In this paper we study multisolitons in a two-dimensional version of the Skyrme modeL 
The model was first considered in 0, but the motivation there is somewhat different from the 
approach taken here. For the present purpose it is important to be clear in what sense our 
model resembles Skyrme's model. In this section we will therefore briefly review a general 
framework for the Skyrme model due to Manton 0] and explain how our model fits into that 
framework. In |^ the usual Skyrme energy functional is interpreted in terms of elasticity 
theory and a static Skyrme field is a map 

7r:5^S (1.1) 

from physical space 5" to the target space S. Both S and S are assumed to be Riemannian 
manifolds with metrics t and r respectively. The energy of a configuration vr is expressed in 
terms of its strain tensor D. To calculate the strain tensor one introduces coordinates on 
S and vr" on E and orthonormal frame fields Sm on S and cr^ on S (1 < i, m < dimS", 1 < 
< dimS). The Jacobian of the map vr is, in orthonormal coordinates, 

and the strain tensor D is defined via 

The energy functional should be a function of the strain tensor but it should not depend 
on the choice of the orthonormal frame s^- The basic invariants under orthogonal trans- 
formations of D, however, are well known: they are the coefficients in the characteristic 
polynomial Xoit) '■= det{D — t id). If S is three-dimensional there are three such invariants, 
namely tiD, |(trD)^ — ^tiD^ and detD. It is explained in that in the usual Skyrme 
model, where S = and S = SU{2), the energy functional is constructed from the first 
two: 



-^S kyrme 



J d^x (iiD + ^(trD)2 - ^trD^^ . (1.4) 



It is now straightforward to construct a Skyrme model in two dimensions in this frame- 
work. We will reserve the term Skyrme model for Skyrme's original model and refer to our 
two-dimensional model as a baby Skyrme model, its soliton solutions baby Skyrmions etc. 
We want to work in flat space, so we choose S = R^. The choice of S is less clear. In order 
to obtain solitons with a topological charge we require, as we shall explain, 7r2(S) = Z, so 
the simplest choice is S = S*^. Thus a baby Skyrme field is a map 

: R2 ^ S^, (1.5) 



2 



where S"^ is the unit 2-sphere in euchdean 3-space with the metric induced by that embedding, 
and we think of as a three component vector (01, 02, (ps) satisfying 0^ = 0^ + 02 + = 1- 
In two dimensions there are only two invariants of the strain tensor, namely 

62 = trD = + 920' 

64 = detD = (dicf) X a2</))^ (1.6) 

where di,i = 1,2, denotes the partial derivative with respect to the cartesian coordinates 
of the vector x G and x is the vector product in three dimensions. 

The simplest choice of the energy functional is E„ = ^ J d?x 62 and leads to the much 
studied two-dimensional a-model whose static solutions are harmonic maps from to 5*^. 
However, the cr-model is not a good analogue of the Skyrme model since, unlike Skyrmions, 
its soliton solutions have an arbitrary scale. The scale invariance is broken by adding the 
Skyrme term | / dPx 64, but the resulting energy functional is still not satisfactory: it can have 
no minima since the energy of any configuration can be lowered by rescaling: (f){x) 4*(^)-i 
where A > 1. Hence it is necessary to include a term in the energy functional which contains 
no derivatives of the field cf) and which is often called a potential. Here we will consider 
the potential /i^(l — n-cj)), where n = (0,0, 1) and /i is a constant with the dimension of 
inverse length. Our potential is analogous to the extra term included in the Skyrme model to 
give the pions a mass (thus l//i may physically be interpreted as the Compton wavelength 
of the mesons in our model). Both in the Skyrme model and in our model the potential 
term reduces the symmetry of the model and is responsible for the soliton solutions being 
exponentially localised in space. In the Skyrme model the size of Skyrmion is slightly smaller 
than the pion's Compton wavelength, reflecting the relative magnitudes of nucleon size and 
the pion's Compton wavelength in nature. To mimic these properties we want to have a 
basic soliton solution whose size is of order l//i. This can be achieved by setting /i^ = 0.1, 
which we do for the rest of this paper. Note, however, that a different choice of fi results 
in a different model. While some of the general features to be discussed in this paper are 
independent of /i, other properties, such as the shape of multisolitons (we give a precise 
definition furhter below), may well depend on it. 

Some readers may want ot bear in mind an alternative physical interpretation the field 
(f): it can also be thought of as the magnetisation vector of a two-dimensional ferromagnetic 
substance 0]. Then the potential term describes the coupling of the magnetisation vector 
to a constant external magnetic field. 

To sum up: the general framework for Skyrme models described above, augmented by a 
potential modelled on the one used in the usual Skyrme model leads to the following energy 
density for the baby Skyrme model: 

e = ^{d,cf>' + 9202) + ^(9i</> X 920)2 + fi\l - n-cf>), (1.7) 
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from which the energy functional is obtained by integration 

E[(f)] = J £xe. (1.8) 

We are only interested in configurations with finite energy, so we define our configurations 
space Q to be the space of all maps (f) : ^ S"^ which tend to the constant field n, called 
the vacuum, at spatial infinity 

lim 4>{x) = n. (1.9) 

As a result, a configuration may be regarded as a map from compactified physical space 
U {oo}, which is homeomorphic to a 2-sphere, to S"^. Thus every configuration may 
be regarded as a representative of a homotopy class in 7r2(S'^) = Z and has a corresponding 
integer degree, which can be calculated from 

deg[0] = — / (fxct)-di(f) X ^20. (1.10) 

Configurations with different degrees cannot be smoothly deformed into each other and hence 
the configuration space Q is not connected. We write Q„ for the component of Q containing 
the configurations of degree n. In this paper we are interested in stationary points and, if 
they exist, minima of in a given sector Qn. A configuration is a stationary point of 
E if the first variation of E under (t){x) ^ <P{x) + e(x) x (t>{x) vanishes for any function 
e : R^ ^ R^ satisfying (f){x)-e{x) = for all x. This requirement leads to the Euler-Lagrange 
equation 

diji = ii^n X 0, (1.11) 

where 

ji = (px di4> + dj(f){dj(f)-(f) X di(f)). (1.12) 



The degree gives a useful lower bound on the potential energy, the Bogomol'nyi bound 

> 47r-|deg[0]|. (1.13) 

This inequality holds already for the cr-model energy functional E^r (for a proof see p ) and 
since E is bounded below by E^j it holds for E as well. Assuming that there are finite energy 
configurations of degree n it follows that the infimum of the restriction E\q^ of the functional 
E to Qn exists, and we call it En. It is not clear, however, whether there is a configuration 
of degree n whose energy is En- If such a configuration exists for n = 1 we call it a 1-soliton 
or a baby Skyrmion, and if it exists for n > 1 and if moreover its energy satisfies 

En < Ek + El for all integers 1 < k,l < n such that k + I = n (1-14) 
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we call it a multisoliton or an n-soliton. If we do not need to specify the degree we simply 
write soliton for an ra-soliton with arbitrary n G N (all solitons of negative degree can be 
obtained from solitons of positive degree via the iso-reflection (0i, 02,03) ^ (0i7— 02,03), 
so we can restrict attention to positive n without loss of generality). We have included the 
condition ( |1.14|) in our definition because we want multisolitons to be stable with respect to 
decay into multisolitons of smaller degree. 

The goal of this paper is to show the existence of multisolitons in the baby Skyrme 
model and to describe their properties. In section 2 we discuss a particular ansatz for 
finding solitons and use it to calculate the field of a baby Skyrmion. In section 3 we set up a 
general framework for proving the inequality (|1.14| ) and in section 4 we complete the proof 
for certain values of k, I and n. Section 5 contains numerical evidence for the existence of 
n-solitons in our model for 1 < n < 6 and a description of their properties. 



2 Hedgehog Fields and Baby Skyrmions 



A powerful method of finding stationary points of the energy functional E ( |1.8| ) exploits the 
invariance of E and Qn under the symmetry group 

G = E2xSO{2),,^xP. (2.1) 

Here E2 is the euclidean group of translations and rotations in two dimensions which acts 
on fields via pull-back. S0{2)i^^ is the subgroup of the three-dimensional rotation group 
acting on S"^ which leaves n fixed. We call its elements iso-rotations to distinguish them 
from rotations in physical space. Elements of 30(2)^^^ can be parametrised by an angle 
X G [0, 27t) and act on via 

(01,02,03) ^ (cos X 01 + sin X 02, - sin X 01 + cos x 02, 03)- (2.2) 

Finally P is a combined reflection in both space and the target space S'^: 

P : (xi,a;2) ^ (xi, -X2) and (0i, 02, 03) (0i, -02, 03) (2.3) 



The vacuum field n is invariant under the symmetry group G and clearly minimises the 
energy in Qq. However, we are interested in stationary points of degree 7^ and the maximal 
subgroups of G under which such a field can be invariant are labelled by a non-zero integer 
n and consist of spatial rotations by some angle a G [0, 27r) and simultaneous iso-rotation 
by —na. Fields invariant under such a group are of the form 

4>{x) = (sin /(r) cos(n6' — x), sin /(r) sm{n9 — x), cos /(r)), (2.4) 
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where (r, 6) are polar coordinates in the cc-plane, and / is function satisfying certain bound- 
ary conditions to be specified below. The angle x is also arbitrary, but fields with different 
X are related by an iso-rotation and therefore degenerate in energy. Hence we concentrate 
on the fields where x = 0- Such fields are the analogue of the hedgehog fields in the Skyrme 
model and we also call them hedgehog fields here. They were also studied in P] for a different 
value of /i. 

The function /, which we call the profile function, has to satisfy 

/(O) = mvr, me Z, (2.5) 



for the field (|2.4|) to be regular at the origin and to satisfy the boundary condition (|1.9| ) we 

set 

;im/(r)=0. (2.6) 

We can assume without loss of generality that m is positive because changing the sign of / 
in (|2.4| ) is equivalent to an iso-rotation by 180°. The restriction E of E to fields of the form 
1) is 

E = 2nJ rdr (^i/'^ + ^^^f^(l + f") + - cos/)) . (2.7) 

It follows from the "principle of symmetric criticality" |^ that a field of the form ( ^.4[ ) is a 
stationary point of the energy functional if / is a stationary point of E, i.e. if / satisfies 
the Euler Lagrange equations first written down in 

^ n^sm'^f\ „ ( n^sinV , n^/' sin/ cos /\ , sin /cos/ 2 ■ c n 
r + / +1 + / r/i sm / = 0. 



(2.8) 



The behaviour of solutions of this equation near the origin and for large r can be deduced 
analytically and was also discussed in 0. The result is that, for small r, 

/(r) ^ mTT + Cr'^, (2.9) 

where C is a constant which depends on m and n. For large r, the equation ( f^.8| ) simplifies 
to the modified Bessel equation 

1 r?2 

r + -/'-(^ + /i')/ = 0. (2.10) 

Solutions of this equation which tend to zero at r = 00 are the modified Bessel functions 
Kn{fir) of order n. Thus a solution / of ( |2.8D is proportional to Kn for large r and we can 
write 

m ~ ^/^„(/ir), (2.11) 
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where c„ is a constant, dependent on n and m, which we will interpret further below. Since 
the modified Bessel functions have the asymptotic behaviour 



we know that the leading term in an asymptotic expansion of / is e~^^'/y/r. The behaviour 
of / at infinity guarantees that the energy distribution of the corresponding field 4> is expo- 
nentially localised and hence that its total energy is finite. 



What is the degree of the field (|2.4| )? A short calculation gives 



deg[(f)] 



n if m odd 
if m even. 



Thus for a given degree n 7^ there are infinitely many solutions of the static field equations 
of the form ('E.4|), one for each odd m. 



Using a shooting method we have solved the ordinary differential equation ( p.8| ) nu- 
merically for a range of values of n and m and computed the energy of the corresponding 
hedgehog field from (|2.7|). We find that, amongst all solutions of the static field equations 



of the form ( |2.4| ) with degree 1 , the field with n = m = 1 has the lowest energy. We write 
for that field in standard orientation (x = 0) and denote the profile function by f^^\ 
We have also numerically tested the stability of the field (j)^^^ against perturbations which 
destroy the rotational symmetry and found it to be stable (we will describe our numerical 
method in more detail in section 5). Thus, in analogy to the Skyrme model, the minimal 
energy configuration amongst all fields of degree one is of the hedgehog form. The profile 
function f^^^ is plotted in figure l.a) and a plot of the energy density e of the field, which is 
given by the expression in round brackets in ( p.7|) , is shown in figure l.b). The actual field 
4>^^'^ is displayed in figure 2. a). 

There is a whole manifold A^i of minima of E\q^ obtained by acting with G on 
All elements of A^i are baby Skyrmions as defined at the end of section 1. TVli is three- 
dimensional, so a baby Skyrmion is characterised by its two-dimensional position vector and 
an angle x specifying its orientation. The value oi E o\i M.i physically represents the mass 
El of a baby Skyrmion. Numerically we find Ei = 1.564-47r. 



3 Existence of Multisoliton Solutions 



The basic ingredient for proving (|1.14|) is a superposition procedure for solitons. In our 
model such a procedure can be found using the stereographic projection 



p: CU {00} 
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This allows us to translate a function 



M : U {00} C U {(X)}, m(cx))=0 (3.2) 
into a configuration 0" via 

0"=poM (3.3) 

and vice versa. If we think of m as a function of z = Xi + 1x2 and z = Xi — 1x2 then its 
degree, which equals the degree of 0", is the number of poles in 2;, counted with multiplicity, 
minus the number of poles in ^, also counted with multiplicity. For more details on the 
formulation of two-dimensional Skyrme models in terms of u instead of (f) we refer the reader 



to |T0|. Here we use it to define: the superposition of the configurations 0" and (p" is the 



configuration (f)^ where w = u + v. Then clearly deg[<;;i)"'] = deg[0"] + deg[</)'']. 

Now let 0" and (f)" be multisolitons of degrees k and /. We want to show that under 
certain circumstances -£'[0"'] < -E[0"] + -E'[0"]. The idea for our proof is taken from an 
unpublished paper by Kugler and Castillejo where the authors claim to to prove the 
corresponding statement for the Skyrme model. However, their proof contains an unjustified 
assumption as we shall see further below. 

Consider the situation where the multisolitons 0" and 0" are well separated. More 
precisely we assume that we can divide into two regions such that u is small in region 
2 and v small in region 1. The smallness of u means that 0" is close to the vacuum n in 
region 2 and hence of the form 



0" = ^1 - (^"-(^"n + yj" ^ n + (p" + C(((p"-c^"), (3.4) 

where c^"-n = 0. Then, if 0" satisfies the Euler-Lagrange equation ( |1.11| ), cp" satisfies the 
linearised equation 

(A - /i2)(p'^ = (3.5) 

in region 2. Similarly, (f)" is close to the vacuum in region 1 and one defines 93^ analogously 
to (^". It satisfies 

(A-/i2)c^" = (3.6) 
in region 1. One checks that 0"" has the following expansion in powers of ^p'" in region 1 

0^" ^0" + e^ X 0" + ^e^ X (e^ X 0"), (3.7) 
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where e" is linear in (p^ but also depends non-linear ly on 0": 



Similarly, in region 2 



= X ((1 + r-n)^" - ir-V^nn)) . (3i 



0"' ^ (/)'' + e" X 0'' + ^e" X (e" X 0''), (3.9) 



where e" depends linearly on 9?" but non- linearly on (p'": 

= X ((1 + 0^.n)c^" - (^'^■cp")/!)) . (3.10) 

Using these formluae we evaluate £'[</)"']: in region 1 we keep terms which involve (f)'^ 
only, terms which are linear in cp^ and those terms quadratic in (p^ which are independent 
of 0". Similarly, in region 2 we keep terms which involve (f>'" only, terms which are linear y?" 
and those terms quadratic in y?" which are independent of 0^. Denoting integration over 
region 1 and 2 simply by the suffix 1 and 2 respectively we find 



E[(f)''] ^ J^d^xe {4)^ + e'' X 0" + ^e'' x (e^ x 0") 

+ j^d^xe (^0^ + e" X 0" + ie" x (e" x 0") 

^ ^rf2xe(0") + j^d\ {^d.if^-d,^^ + ^/iV"-^'") 

+ ^d2^e(0^) + Q^.v^^-^.v?^ + ^/iV^-¥'^) 

+ ^rf2xj7-9ie'' + /x2e-.nx 0" + ^rf'xj7-9ie" + /i2e"-nx 0^ (3.11) 



where jf" and jl denote the current jj (|1.12|) evaluated on the fields 0" and 4>' respectively. 



The integrals in the third line represent the contribution to energy of the soliton 0" from 
region 1 and the leading contribution from region 2. The integrals in the third line are the 
contribution to the energy of cf)" from region 2 and the leading contribution from region 1. 
The formulae for subleading contributions are complicated, but it is clear that the sum of 
all terms which only involve either 0" or (f)^ is just -E[0"] + -E[0^]. The cross terms are more 
interesting. Integrating by parts in the last line of the equation above, and using the fact 



that both 0" and cf)^ satisfy the Euler-Lagrange equations ( |1.11| ), we are only left with a 
boundary term 

E[r] ^ Em + Em + jU--e- - jrendS,. (3.12) 

Here F is a curve without self-intersections separating the region 1 from the region 2 and 
dSi = Eij'jjdt for any parametrisation of 7(t) of F for which region 1 is on the left and region 
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2 on the right {Sij is the antisymmetric tensor in two dimensions normahsed so that £12 = 1). 
We now assume further that F hes in a region where both (f)'^ and (f}"" are close to the vacuum 
and keep only terms which are linear in both cp^ and (p"". Then, interpreting the difference 
E[(f)'^] — E[(f)'^] — E[(f)'"] as the potential describing the interaction of the solitons 0" and (j)'" 
we find that the leading term V of that potential is given by the simple formula 

V = J^iif'^-d^ip'' - cp''-d,ip^)dS,. (3.13) 



To prove ( |1.14| ) one needs to show that we can always arrange for V to be negative. 



Suppose that this is not already the case, and consider the case where V is positive. Then it 
can be made negative by an iso-rotation by 180° of either 0" or (/>'" (which does not change 
the individual energies of 0" and 4>^). However, the proof is incomplete unless we can rule 
out that V is zero. This possibility was not considered in [^. In fact we shall see that 
this does indeed happen for all relative iso-orientations in a slightly modified version of our 
model. Thus we can only use the result ( |3.13|) for constructing new multisolitons out of 
solitons for which the asymptotic field if is known. 

4 New Multisolitons from Old 

Although we already have enough information about the asymptotic field of a baby Skyrmion 
to prove the existence of a 2-soliton in our model, it is useful to interpret the asymptotic 
field further before proceeding with the proof. 

For large r the profile function f^^^ approaches exponentially fast and we can therefore 
approximate sin/*^^^ ~ f^^'^ and cos/*-^-* ~ 1. Using the asymptotic expression ( 2.11|) we 
write the asymptotic form c^*-^) of (p^^^ as 

= ^(K.ifir) cos(^ - x), i^i(/ir) sm{e - x), 0), (4.1) 
where we have written p for ci. Alternatively, introducing the orthogonal vectors 

Pi = p(cosx,sinx) P2 = P(-sinx,cosx) (4.2) 
and X = x/r we can write 

^«(^) = ^p^-xK,{fir) = -^p^.VKoifir) a = 1,2. (4.3) 

ZTT ZTT 

Now, given that the Green function of the static Klein-Gordon equation in two dimensions 
is Ko(/ir): 

(A - i2'^)Ko{i2r) = -27r6^^\x), (4.4) 
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it follows that 



(A - ^?)^'i\x) = p^-V6^'\x) a = 1,2. (4.5) 

Thus the asymptotic field cp^^^ may be thought of as produced by a pair of orthogonal dipoles, 
one for each of the components (pi'^ and (p2'\ in a linear field theory, namely Klein-Gordon 
theory. The strength of the dipole can be calculated from the asymptotic form of f^^\ We 
find 

p = 24.16. (4.6) 

Consider now the set up of section 3, with 0" a baby Skyrmion centred at the origin 
and iso-rotated relative to the standard hedgehog field (p^^^ by Xi and (f)^ a second baby 
Skyrmion iso-rotated relative to the standard hedgehog by an angle X2 and centred at R, 
where R .= \R\ << 1/yU. The asymptotic field of the first baby Skyrmion is 

= -^da-VK,{fir) a = 1,2. (4.7) 

where 

di = p(cosxi,sinxi) 0^2 = p(- sinxi, cosxi), (4.8) 
and the asymptotic field of the second is 

= -^Pa-^KM^ -R\) a = 1, 2, (4.9) 

where 

Pi = p(cosx2,sinx2) P2 = P(-sinx2,cosx2)- (4.10) 

Then, using that, in region 1, 

{/^- 1X^)^1 = da-V5^^\x) and (A - /i2)v?^ = a = 1,2 (4.11) 

and converting the line integral ( p. 131) into an area integral over region 1 we find the potential 
for the interaction of two well separated baby Skyrmions 



y: ^(d,-v)(p,-v)iro(/ii?) 

a=i,2 2vr 

— cosipAKo^fiR) 

TT 
2 2 

^^cosijKoifiR), (4.12) 

TT 
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where ip = Xi ~ X2 describes the relative iso-orientation. In the last step we have used ( [4 .41 ) 



and have omitted the ^'^^^-function term because we are only interested in large separations 
R > Thus we see that the potential for the interaction of two well separated baby 

Skyrmions is the same as that calculated in a linear field theory, namely Klein Gordon 
theory for a pair of scalar fields, for the interaction of two pairs of orthogonal dipoles. The 
important point is that baby Skyrmions not only act as sources of dipole fields but also react 
to an external field like a pair of orthogonal dipoles. This result deserves some comments. 
Firstly, it is instructive to compare it to a similar result in the Skyrme model. There a 
certain non-linear superposition procedure, the product ansatz, also leads to an interaction 
potential which has a "linear" interpretation: it describes the interaction between two triplets 
of mutually orthogonal scalar dipoles in a linear field theory for the pion fields, see and 
also 0. 

Secondly, we are now in a position to explain the caveat at the end of the previous 
section concerning the use of ( p.l3| ) for proving ( p..l4| ). In |]I0| a similar model to ours was 
studied with the potential term — <p-n) replaced by (1 — <p-n)'^. In that model (where 
the mesons are massless) the solitons of degree 1 also have the hedgehog form ( p.4| ) (with 
a particularly simple profile function) and their asymptotic field can still be interpreted in 
terms of a pair of dipoles. The Green function of the linearised theory, however, is InR. 
Thus a calculation analogous to ours yields a potential proportional to cos ipA In R which is 
identically zero for all values of ip. Numerically, two solitons of degree one are found to repel 
each other so that there are no 2-solitons in that model. The example shows that the "linear 
forces" between solitons, calculated via ( p.l3|) , may vanish for all relative iso-orientations. 
In that situation the inter-soliton forces are entirely due to non-linear effects and cannot be 
calculated with the methods of the previous section. 

Returning to the formula ( [4.12| ) we see that Vu is negative if one baby Skyrmion is iso- 
rotated by 180° relative to the other, i.e. iiip = it in the above expression. Thus we conclude 
that E2 < 2Ei. In fact one finds that already the minimum of the energy amongst hedgehog 
fields of degree 2 is less than 2Ei. We write (f)^'^^ for the hedgehog field ( p.4| ) with = 2 in 
standard orientation (x = 0) whose profile function satisfies (|2.6|),(p78D and ( p.5| ) with m = 1. 
For its energy, or mass, we find E2 = 2.936 ■47r. We have again checked the stability of (f)^'^^ 
against more general perturbations numerically and conclude that it minimises the energy 
amongst all fields of degree 2. Its profile function is plotted in figure l.a) and its energy 
density in figure l.b). Note that the maximum of the energy density is not at the origin but 
at r ^ 1.8. This is again reminiscent of the Skyrme model, where the energy of the static 
solution of degree 2 is concentrated in a toroidal region |jll|. The field </)*-^-' is shown in figure 
2.b). 

Before we can use the 2-soliton as an input for the construction of higher multisolitons 
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we need to understand the asymptotic form tp^"^^ of the field (Z)*-^-*. Since the 2-sohton is of 
the hedgehog form this is not difficuh. The asymptotic field is now 

2 

^(2) ^ |^(K2(/ir) cos(20 - x), K2iixr) sin(20 - x), 0) (4.13) 

and can be expressed in terms of the quadrupole moments 

qi = g(cos X, sin x) 92 = sin x, cos x) (4.14) 

and the second order differential operator D = {df — 9|, 2did2) via 

^^^^ = ^q,-DKoifir) a = 1,2. (4.15) 

Then, using again ( |4.4| ), it follows that 

iA-^^')^^^^ = -q,.D5^'\x). (4.16) 

Thus we may think of the asymptotic field (f^"^^ as being due to a pair of orthogonal 
quadrupoles (in two dimensions all multipoles have two real components). For the strength 
of the quadrupole we find 

q = 53.6. (4.17) 

Now consider the superposition of a baby Skyrmion at the origin and iso-rotated relative 
to the standard hedgehog by xi and a 2-soliton centred at -R, i? >> 1/yU, and iso-rotated 
relative to the standard hedgehog 0'-^^ by X2- Thus the asymptotic field and the dipole 
moments of the baby Skyrmion are as in ( |4.7| ) and ( [4. 81 ) and the asymptotic field of the 
second is 

= ^qa-DKoifi\x -R\) a = 1, 2, (4.18) 

where 

^1 = ?(cosx2,sinx2) 92 = ?(-sinx2,cosx2)- (4.19) 



Then, inserting these expressions into the general formula (|3.13|) and using ( [4.1 Ij ) we find 



the potential V12 describing the interaction between a baby Skyrmion and a 2-soliton : 

a=i,2 2vr 

= —(cos ipdi + sin ip ^2) AXo(y^i?) 
Ztt 

3 

cos(V'-?9)i^i(/ii?), (4.20) 



27r 
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where {R, "i?) are polar coordinates for the relative position vector R and ip is defined as 
before. Thus the interaction between a baby Skyrmion and a 2-soliton depends both on 
their relative iso-orientation and on their relative position in physical space. By choosing 
ip — -t} = we can again arrange for the potential to be negative, so we expect there to be 
a 3-soliton solution in our model. Before describing that solution in the next section we 
calculate the potential V22 for the interaction between two 2-solitons. 

Thus we look at the superposition of a 2-soliton at the origin and iso-rotated relative 
to the standard hedgehog field by Xi and a second 2-soliton centred at R and iso-rotated 
relative to the standard hedgehog field by X2- For the former the quadrupole moments are 

ei = g(cosxi,sinxi) 62 = g(- sin xi, cos Xi) (4.21) 



and for the latter the quadrupole moments are and ^3 defined in ( [4.191 ). Thus, the 
interaction potential is 

a=i,2 2vr 



^2 

-—cos'^A'^KoifiR) 



TT 
2 4 

= -ULcosipKoifiR). (4.22) 

TT 

Thus V22 has the opposite sign from Vn, but the same functional dependence on the separa- 
tion R and the relative iso-orientation ip. This deserves a comment, as it may seem surprising 
that the dipole-dipole potential Vu does not necessarily dominate over the quadrupole- 
quadrupole potential V22 at large R. However, it is characteristic of linear field theories with 
an exponentially decaying Green function that the leading term in a field produced by an 
rz-pole is independent of n. This observation will be important for us later, when we study 
multisolitons whose fields are not of the simple hedgehog form; it means that all multipoles 
in the expansion of the asymptotic field are potentially equally important when studying the 
interaction of such multisolitons. 

For now the most important feature of V22 is that it is negative when ip = 0. Thus we 
conclude E4 < E2 + E2. However, without further insights into the properties of the 3-soliton 
we cannot say anything about the relative size of i?4 and Ei + E^. 



5 Numerical Results for Higher Multisolitons 



There are a number of ways to search numerically for stationary points of energy functionals 
and to check whether a given stationary point is a local minimum. It is very hard, however, 
to ascertain whether local minima found in this way are global minima. We are similarly not 
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able to prove that the configurations to be described in this section are multisohtons in the 
strict sense of the definition at the end of section 1. Instead we will consider a numerically 
found configuration to be a multisoliton if it satisfies certain numerical checks for a local 
minimum and the inequality (|1.14|) . 

Here we look for stationary points of the energy functional E by solving a suitable time- 



dependent equation which reduces to (|1 . 1 1| ) in the static limit. The time evolution according 
to that equation should stop at stationary points of E. The equation we use is the Lorentz 
covariant equation 

d^icf) X da(f> - 9^0(9/30-0 X da(f>) = 12"^ (f) X n, (5.1) 

where is now a function of = {t,x) and the indices a, (3 = 0,1,2 are raised and 
lowered with the Lorentzian metric diag(l, —1, —1), with an added friction term to absorb 
the kinetic energy. To solve the resulting equation numerically we use a finite difference 
scheme to evaluate the space derivatives and integrate the time evolution using the 4th 
order Runge-Kutta method. We use a square grid of 200 x 200 points and set the time 
increment dt to half the length of a lattice site. Most of our simulations are performed on a 
grid extending in both the and direction from —20 to 20. 

As the initial configuration we take a particular stationary point of E, namely the hedge- 
hog field ( p.4|) of degree n with the profile function satisfying the boundary condition ( p.5|) 
for 171 = 1 and the ordinary differential equation ( p.8|) . We then add a small perturbation 
which breaks the symmetry of the hedgehog. We have already noted that the hedgehog fields 
for n = 1 and n = 2 are stable against such perturbations. We have also investigated the 
cases 3 < n < 6: in those cases the hedgehog fields are unstable and the time evolution ends 
at less symmetric configurations with lower energy. These final configurations are stable 
with respect to further perturbations, and their energies satisfy the inequality ( p..l4|) so we 
take them to be multisohtons. 

From the energies E.^ of those multisohtons we calculate the "ionisation energies" A^i 
defined via 



^kl 



Ek + Ei- En, where 1 < A;, / < n, k + I = n. (5.2) 



In table 1 we summarise our results. All the energy values listed in the table are obtained by 
integrating the solitons' energy density over the grid. The resulting values for Ei and E2 are 
1% smaller than the more accurate values given earlier in the paper, which were obtained by 
solving the ordinary differential equation ( ^^ and integrating (|2.7| ). For the calculation of 
the ionisation energies, however, it is important that all energies are calculated in the same 
way. 
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n 


E„/4:TT 


k + l 


AwMtt 


1 


1.549 






2 


2.907 


1 + 1 


0.191 


3 


4.379 


2 + 1 


0.077 


4 




3 + 1 


1 28 






2 + 2 


0.014 


5 


7 282 


4 + 1 


068 






3 + 2 


005 


6 


8.693 


5 + 1 


0.138 






4 + 2 


0.015 






3 + 3 


0.066 



Table 1 

Multisoliton energies and ionisation energies 



The results for n > 3 deserve a more detailed discussion. For n = 3 our numerical 
procedure leads to the configuration displayed in figure 3. a). A plot of the energy density 
e ( |1.7| ) as a function of position is shown in figure 3.b). Unlike the solutions discussed so 
far the energy of the 3-soliton is not rotationally symmetric. Instead the configuration is 
like a linear molecule made up of three (distorted) baby Skyrmions aligned so that any two 
neighbours are in the most attractive relative orientation. 

In principle it is still possible to analyse the asymptotic field of the 3-soliton in terms of 
multipole moments. In practice the absence of a continuous symmetry makes the analysis 
hard and we know from the previous section that we may have to consider many multipole 
moments. We have therefore made no attempt at deriving the potential for the interaction 
between 3-solitons and other multisolitons in our model and rely on numerical evidence for 
showing the existence of higher multisolitons. 
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The field and the energy density for the 4-soliton are plotted in figure 4. The 4-soliton 
is again like a linear molecule, but this time made up of two 2-solitons. The plot of the field 
shows that the two 2-solitons have the same iso-orientation, as expected from our discussion 
of the potential ¥22- The picture suggests, and table 1 confirms, that it costs very little 
energy to break the 4-soliton into 2-solitons. 

The field for the 5-soliton is the least symmetric of all the multisolitons we have studied. 
The field and its energy density are shown in figure 5. The 5-soliton consists of an almost 
undistorted 2-soliton and 3-soliton close together. Table 1 shows that the binding between 
those constituents is very weak. To check whether the binding is in fact a boundary effect 
we have therefore repeated the simulation on a grid of 250 x 250 points extending in and 

from —25 to 25. The result is identical to that of the first simulation. 

Finally the 6-soliton is made up of three 2-solitons centred at the vertices of an equilateral 
triangle. The plot of the field in figure 6. a) reveals that the 2-solitons all have the same iso- 
orientation, so that the interaction energy between any two 2-solitons is minimised. Although 
the binding is again quite weak, the distortion of the individual 2-solitons as a result of their 
interaction is clearly visible in the plot of the energy density in figure 6.b). 

The pictures of the energy density for n > 4 suggests that the 2-soliton (and to a lesser 
extent the 3-soliton) serves as a basic building block for higher multisolitons. Table 1 provides 
further evidence for this observation: the largest ionisation energy in table 1 is An, showing 
that the 2-soliton is most strongly bound, and the ionisation energy A^/ for n = k + I > 2 
is least if A; = 2 or / = 2, showing that it is easiest to break up multisolitons in a way that 
produces at least one 2-soliton. 



6 Conclusions 



We have studied multisoliton solutions in a non-linear field theory which may be interpreted 
physically as a two-dimensional version of the Skyrme model for nuclear physics. Our main 
analytical result is the expression ( p.l3| ) for the leading term in the potential describing the 



interaction of two well separated solitons in terms of the asymptotic fields of the solitons. 
In those cases where the asymptotic field of the solitons is known the leading term could 
be written down explicitly in terms of the multipoles associated with the solitons and the 
Green function of the linearised theory. The resulting formula can then be used to prove the 
existence of multisolitons of higher degree. 

Using numerical methods we could explicitly display multisolitons of degree 1 < n < 6. 
It turns out that in our model solitons of degree 1 and 2 are invariant under a SO (2) action, 
but that higher multisolitons are only invariant under the action of finite groups. This is 
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rather reminiscent of the Skyrme model where sohtons of degree 1 and 2 have continuous 



symmetries |]TT|, but higher muhisohtons have only discrete symmetries, see |T2[ and also 

Hi- 

The relationship between solitons and the associated multipoles is intriguing and deserves 
further study. It is quite possible that one can prove general statements about the multipole 
expansion of the asymptotic field of an n-soliton for arbitrary n which would allow one to 
complete the proof of the inequality ( |1.14|) for general /c, / and n. This would be of wider 
interest since the method described in section 3 applies to any field theory with soliton 
solutions provided there is some sort of superposition procedure for well-separated solitons 
and the linearisation of the theory is of a suitable form. 

Finally it would be interesting to see to what extent the multipole description can be 
used to understand the dynamical properties of the solitons in our model. The form of the 
potential Vu shows that the forces between two well separated baby Skyrmions, like those 
between Skyrmions, depend both on the relative separation and the relative orientation. This 
suggests that the interactive dynamics of two baby Skyrmions is more complicated than that 
of other topological solitons in two dimensions, such as CP^ lumps [0] or the solitons studied 
in |10|, but also possibly more relevant to Skyrmion dynamics in three dimensions. We are 



presently investigating this point and will report on it elsewhere |T5 
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Figure Captions 



Figure 1 

a) Profile functions for tlie baby Skyrmion (n = 1) and the 2-soliton {n = 2). 

b) Energy densities e ( |1 . 7| ) as a function of r for the baby Skyrmion (n = 1) and the 2-sohton 
{n = 2). 

Figure 2 

a) Plot of the field (()^^\ At every lattice site in physical space we plot an arrow whose 
direction and magnitude is that of {4>i'\4'2^^) (we identify the axis in the target space S'^ 
with those in physical space). At the base of the arrow we put a '+' if 03^'' is positive and a 
'x' if 03^^ is negative. Thus the vacuum is represented simply by a '+'. The labels x and y 
refer to the first and second component of the vector x. 

b) Plot of the field (f)^'^^ using the same conventions as in a). 

Figure 3 

a) Field of the 3-soliton; conventions as for figure 2. a). 

b) Energy density of the 3-soliton in the range —10<x,y< 10. 

Figure 4 

a) Field of the 4-soliton; conventions as for figure 2. a). 

b) Energy density of the 4-soliton in the range —10 < x,y < 10. 

Figure 5 

a) Field of the 5-soliton; conventions as for figure 2. a). 

b) Energy density of the 5-soliton in the range —10<x,y< 10. 

Figure 6 

a) Field of the 6-soliton; conventions as for figure 2. a). 

b) Energy density of the 6-soliton in the range —10 < x,y < 10. 
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